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Abstract. The inverse seesaw mechanism allows the neutrino masses to be generated by 
new physics at an experimentally accessible scale, even with 0(1) Yukawa couplings. In 
the inverse seesaw scenario, the smallness of neutrino masses is linked to the smallness of 



^S^ a lepton number violating parameter. This parameter may arise radiatively. In this paper, 

we study the cosmological implications of two contrasting radiative inverse seesaw models, 

1-H one due to Ma and the other to Law and McDonald. The former features spontaneous, the 

^ latter explicit lepton number violation. First, we examine the effect of the lepton-number 

\^ violating interactions introduced in these models on the baryon asymmetry of the universe. 

We investigate under what conditions a pre-existing baryon asymmetry does not get washed 

out. While both models allow a baryon asymmetry to survive only once the temperature 

^r has dropped below the mass of their heaviest fields, the Ma model can create the baryon 

~^ asymmetry through resonant leptogenesis. Then we investigate the viability of the dark 

T— I matter candidates arising within these models, and explore the prospects for direct detection. 

^ We find that the Law/McDonald model allows a simple dark matter scenario similar to the 

Higgs portal, while in the Ma model the simplest cold dark matter scenario would tend to 

overdose the universe. 
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1 Introduction 

The experimental evidence for physics beyond the standard model (SM) includes the baryon 
asymmetry of the universe (BAU), the existence of dark matter (DM), and the non-zero 
neutrino masses. It is possible that the mechanisms responsible for these observations are 
intrinsically linked. For example, the observed neutrino masses may arise from new interac- 
tions which at low energies generate the Weinberg operator [1], Cw = (1/A) ^^1^1^, where 
<1> is the SM Higgs, II are the SM lepton doublets and A is the scale of new physics. Cw is 
the lowest-order (non-renormalisable) operator containing only SM fields which can give rise 
to neutrino masses. Moreover, it violates the lepton-number (L) symmetry of the SM, thus 
suggesting that the interactions which give rise to it may have also generated the BAU via 
leptogenesis at early times [2]. It is also plausible that the relic abundance of some of the 
degrees of freedom mediating Cw constitutes the DM of the universe today [3-9]. 

Another possibility is that the BAU was generated by a mechanism unrelated to the 
physics of neutrino masses. Well-motivated scenarios include the Affleck-Dine mechanism 



in supersymmetric models [10, 11], and baryogenesis via a first-order phase transition [12]. 
Clearly, whatever mechanism of baryogenesis operated in the early universe, the existence 
of a baryon asymmetry today has to be consistent with the interactions which give rise 
to the neutrino masses. If the latter violate L, then they could pose a challenge for the 
BAU: L-violating interactions operating efficiently in the early universe, together with the 
rapid electroweak (EW) sphaleron processes [13, 14], which violate the baryon-plus-lepton 
number of the SM, can erase any baryon-number (B) asymmetry created before the EW 
phase transition. It is then important to understand what conditions the survival of a baryon 
asymmetry imposes on neutrino mass generation mechanisms. 

The three minimal, tree level, UV completions of the Weinberg operator, the type-I [IS- 
IS], type-II [19-24], and type-Ill [25] seesaw mechanisms, have been studied extensively for 
their phenomenological consequences. In these scenarios, the scale of new physics ranges 
from ~100 GeV for Yukawa couplings of 0{W~^), all the way up to ~ 10^^ GeV for Yukawa 
couplings of 0(1). The small Yukawa couplings in the former case, and the very high energy 
scale in the latter case make these models difficult to test experimentally (although for the 
type-II and III seesaws gauge interactions allow the lower mass range to be tested at the 
LHC). 

On the other hand, the inverse seesaw (ISS) mechanism [26, 27] can generate the ob- 
served neutrino masses via new physics at an experimentally accessible scale, and with 0(1) 
Yukawa couplings. In the ISS, the smallness of neutrino masses is traced to the smallness of 
an L-violating parameter, which we shall denote by //. This small lepton-number violation 
naturally suppresses the washout of a pre-existing baryon asymmetry [28]. To explain why /x 
is small, it is natural to suppose that it is generated radiatively [29-33]. (Other possibilities of 
explaining why // is a small parameter are extra dimensions [34, 35] , or a type-II seesaw-like 
mechanism [36, 37].) 

Radiative ISS models introduce new fields and lepton-number violating interactions, 
which could potentially wash out the BAU. In this paper, we investigate the constraints on 
the parameters of the ISS implied by the BAU, or conversely, the conditions on baryoge- 
nesis assuming that neutrino masses are generated via the ISS. We consider two radiative 
ISS models: one in which lepton number is explicitly broken [32], and one in which it is 
spontaneously broken [29]. These models also feature accidental discrete symmetries which 
stabilise exotic particles. The latter can play the role of DM. We investigate the parameter 
space, derive conditions which ensure that the relic abundance of the new stable particles 
does not overdose the universe, and identify more specifically the regions which give rise to 
the observed DM abundance. Similar issues in a different implementation of the ISS have 
been studied in ref. [38]. 

The order of the paper is as follows. In section 2 we introduce the ISS, and review 
the natural suppression of washout in the ISS. In section 3 we introduce the Law/McDonald 
radiative ISS model [32], analyse washout with the additional fields present, and examine its 
DM candidates. In section 4 we repeat the analysis for the Ma radiative ISS model [29]. 

2 The inverse seesaw mechanism 

2.1 The generic mass matrix 

The ISS [26, 27] is generically realised with the additional terms in the Lagrangian 

C D X^e'^'h^^lNR + MrN^Sl + l^SLTiSL) + J/i(iV^(iVR) + H.c, (2.1) 




where $ ~ (1,2,1/2) is the usual SM Higgs doublet, II ~ (1,2,-1/2) is the SM lepton 
doublet, and Nr and Sl are singlets under the SM gauge group. We leave the details of 
additional symmetries and symmetry breaking which results in the above terms [39] unspec- 
ified for now (a field redefinition ensures the absence of the lL^^{SLy term). After $ gains 
a vacuum expectation value (VEV) ($) = (0 , f^/\/2)"^ and breaks the EW symmetry 
SU{2)l ® U{1)y — )• C^(1)em, the neutral-fermion mass matrix is (assuming one generation) 

-, ,0 ruD 

C = -[l7ZN-^SL)\mD A Mr\\ Nr I , (2.2) 

Q Mr ^l 

where rriD = ^uVw/v^- Diagonalization of the mass matrix for n,fl <^ itij^i^Mr leads to a 
light neutrino mass 

2 
777 

-. « ^/^ . (2.3) 

-ft 

The heavy Majorana neutrinos A'^i and A''2, superpositions of mostly (SlY and Nr, have 

mass eigenvalues 

Mn^Mr±'^. (2.4) 

Due to the small mass splitting, they form a pseudo-Dirac pair. The mixing between the 
light and heavy states is approximately 9 ~ ttio/Mr. In the limit ;U — )• the light neutrinos 
become massless at tree level, although loop corrections give finite contributions to their 
masses unless both /i,/i — )• [40]. In this limit one also recovers a global lepton number 
symmetry in the Lagrangian and hence small choices of fi, and jl are technically natural. 
This is the small parameter generated radiatively in the ISS models we return to below. 

This discussion can be generalised to a realistic three generation case. The entries in 
the mass matrix (2.2), are promoted to 3 x 3 matrices, and the light neutrino mass matrix is 
then to first order, 

M^ = MdM^^^x{mI)-^mI. (2.5) 

Various studies on the collider phenomenology [41, 42], flavour violating processes [43-46], 
and neutrinoless double beta decay [47, 48] of the ISS have appeared in the literature. 

2.2 ISS and baryogenesis 

Before the EW phase transition, the non-perturbative sphaleron processes rapidly wash out 
any existing B + L asymmetry, while they preserve B — L. However, if there exist other rapid 
L violating processes, then both B and L will be driven to zero. In the ISS, the SM lepton 
number is violated by the parameters /i and Jl. Provided that these parameters are small 
enough, the B asymmetry is not washed out. 

Let us examine the possible washout of a -B asymmetry created at a temperature above 
the heavy neutrino mass scale, Tbg > -^Af, where Tbg is the temperature of baryogenesis. 
In this case, the inverse decays to on-shell heavy Majorana neutrinos Z^, + $ — )■ iVj are not 
Boltzmann suppressed. Such inverse decays followed by decays Ni ^ Il + <I>* can change the 
lepton number by two units. This washout process is shown in figure 1. 

It was found in ref. [28] that in the ISS such washout processes can be naturally 
suppressed, and we review the findings here. For simplicity we take from now on /i = 0, as 
this parameter does not enter the neutrino mass eigenvalues to first order, and focus on /x 
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Figure 1. AL = 2 process which can wash out the baryon asymmetry at high temperatures. Both 
heavy mass eigenstates A^i and N2 contribute, but destructive interference for small mass splittings 
/i ^ M^v crucially limits the washout rate. 



instead. The more in-depth discussion of [28] allows both parameters to be non-zero with no 
major differences to the discussion here. 

For typical mass choices the heavy neutrinos can decay into a component of the Higgs 
doublet and a lepton N — )■ (j) + II- For heavy neutrino masses well above the Higgs mass, 
Mj^ » M/i, the decay rate is 

\2 

Yd = ^Mn, (2.6) 

ovr 

and for temperatures T > M]\[ one includes a Lorentz factor (Mn/T). In the type I seesaw, 
the effectiveness of washout is well captured by the washout parameter [49] 



K 



H 



(2.7) 



T=Mn 



where H ~ g^''^T^/Mpi is the expansion rate of the universe, and g is the number of rela- 
tivistic degrees of freedom. Typically washout processes are ineffective for K < 1, while any 
net B is erased for K ^ 1. In the intermediate regime if ~ 1 — 100 a more careful analysis 
is required. For our purposes however, we adopt the approximation that K < 1 is safe from 
washout. 

In the ISS, however, the small mass splitting between the heavy pseudo-Dirac neutrinos 
leads to destructive interference in the scattering process depicted in figure 1. This leads to 
a suppressed washout rate [28] 



K, 



eff 



KS^ 



/^ 



TdH 



(2. 



T=M. 



where 6 = fi/Tu. As expected, in the limit of lepton number conservation, ^u — t- 0, the 
heavy neutrinos form a Dirac pair, there is complete destructive interference in the AL = 2 
scattering process, and washout of the baryon asymmetry does not occur. Requiring i^cft ^ 1 
to avoid washout translates into a bound: 



A* < Ai. 



-TW "'^^^■ 



(2.9) 



For Ajy = 0.1, AIn = 1 TeV, light neutrino masses mj^ = 0.1 eV are explained by ^ = 330 eV 
(cf. eq. (2.3)). So the light neutrino masses can be explained while avoiding washout of the 



BAU. Alternatively, one can express the above bound in terms of the mixing angle: 

To satisfy current experimental constraints (e.g. from universality of the weak interaction 
and rare leptonic decays) requires roughly, for Mjv above the EW scale, a mixing angle 
Q < O(10~^) [44, 50]. (The usually stringent constraint on the active-sterile mixing angle 
from neutrinoless double beta decay [51] is greatly weakened due to the small lepton number 
violation in the ISS [47, 48].) Future experimental tests at the Q ~ 10~^ level can therefore 
help test the parameter space corresponding to washout avoidance [52]. 

For baryogenesis temperatures below the heavy neutrino scale, Tbg < -W^Af, the number 
density of heavy states is Boltzmann suppressed. This means the thermally averaged rate of 
the net washout process /l + <I> — )• A^j — ?■ //, + <1>* is suppressed by both a Boltzmann factor 
and the small ^u. In this regime one is even safer from washout. 

What about baryogenesis itself in the ISS? The obvious candidate is resonant leptoge- 
nesis [53]. The resonance itself is due to an almost degenerate pair of heavy neutrinos and 
exactly such pairs appear in the ISS. Detailed studies have shown resonant leptogenesis is 
indeed possible in the ISS framework [54], and it is even possible to link this with the DM 
abundance [55]. 

To explain the smallness of [i we now turn to models where it is generated radiatively. 
Because of the smallness of /i, washout of the BAU is suppressed in ISS models. However, to 
generate a finite value of [i radiatively requires the addition of new fields and L-violating in- 
teractions. We will see the effects of these new fields and interactions on washout of the BAU, 
and study the DM candidates in the radiative ISS models. We study the Law/McDonald 
[32] and Ma radiative ISS models in turn [29]. 

3 Law/McDonald radiative inverse seesaw 

3.1 Review of the model 

The Law/McDonald radiative ISS model [32] employs explicit violation of {B — L), so that 
processes other than scattering through the heavy neutrinos crucially contribute to washout 
of the BAU. The model introduces an additional abelian gauge symmetry to the standard 
model gauge group 

G = GsM^U{l)d. (3.1) 

All the SM particles are neutral under this new charge. The following exotic fermions are 
introduced (with convention Qem = -^3 + ^): 

^R,L~(1,3,-1)(0) Ar^,,L~(l,l,0)(l). (3.2) 

Along with the SM Higgs $ ~ (1, 2, l/2)(0), the following scalars are introduced 

e~(l,3,-l)(-l), r?~(l,l,0)(l), x~(l,l,0)(2). (3.3) 

The following Yukawa interactions and bare Dirac mass terms are allowed 

-C D +yETL^R + hE^XNR + h'^E^^NL + h^N^N'^x + h'J^N^x .3 ^. 

+MeElEr + MnNlNr + H.c. ^ ■ ' 



'X 



2/0 . 



The VEV pattern is chosen to be 

(/) = t;^/%/2 = 174 GeV, (e"), (t?) = 0, {x) = n 
The tree level mass matrix for the neutral fermions is then 

VEVw/y/^ \ ( {vl) 



C 



1 



uZ [ElY El) I yEVy,/V2 





(3.5) 



(3.6) 





Me y} / w^L) 

Note Nl^pi are heavy Majorana fermions decoupled from the above matrix, and we have not 
yet taken into account the radiative corrections which will generate the terms /x(£'^)^-E^ and 
^{E^YEI after symmetry breaking. The scalar potential is given by, 

Vs = {^^^) [4 + X^{^^^) + A^^(ete) + A^^r/tr?) + \4,^{xh) 

+ (e^O [4 + hii^O + kviv^v) + kxix^x) 

+ (^^^) [m' + \{V^V) + >^vxixh)] (3.7) 

+ ixh) \fil + x^ixh) 



+ 



hct>r, C^^V + 2^™ VVX^ + H.c. 
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Figure 2. Radiative mass generation in the Law/McDonald model. 

The SM lepton number is violated explicitly by the combination of eq. (3.7) and eq. (3.4), 
but a Majorana mass term for the El and E^ fields is forbidden due to hypercharge. After 
X and <1> gain VEVs, a radiative mass term fi{ElYEl is generated, as shown in figure 2. A 
mass term fl{EjfYEl^ is also generated radiatively. The radiative mass is approximately 



/^ 



1 HKH^yvivj 
167r2 120 Afgg 



^'?X' 



(3.8) 



where Ajss ~ Mn, Mr,, M^ is the scale of the new physics, usually taken to be ~ ©(TeV), and 
one may choose firix ~ ©(TeV). This radiatively generated /i is exactly the small parameter 
required in the mass matrix for the ISS. The light neutrinos masses are then 

(yppx^^\ (v^\ 

\ 10-5 ) \Me) V A 



^ /I TeV\^ 



rrip 



ISS 



V 



f^vx 
ITeV 



X 0.1 eV 



(3. 



3.2 Constraints from BAU washout 

Previously we have seen that in the ISS framework, the washout of the BAU is natu- 
raUy suppressed due to the smaUness of the lepton-number violating parameter //. In the 
Law/McDonald radiative ISS model, // is in fact generated only after EW symmetry breaking 
and cannot itself cause washout. However, the SM lepton number is explicitly violated by 
the couplings which eventually generate fi and give rise to the radiative ISS. The interactions 
of eqs. (3.7) and (3.4) can then lead to washout. Here, we investigate under what conditions 
the BAU is not washed out. 

The {B—L) of the SM remains a good symmetry if any of the following sets of parameters 
vanishes 

{ve}, 1/^^}' {h<t>v}^ ih^ ^^i' "L^x' ^xi • (^-10) 

(Note that Me has to be large since Er^l are charged under the SM gauge group.) The BAU 
can change only if interactions which involve all of the above couplings are efficient. We show 
such a sequence of interactions in figure 3. In order to preserve the BAU, we require that the 
interactions mediated by at least one of these sets of couplings are out of equilibrium in the 
period between baryogenesis and the EW phase transition. We discern two cases, depending 
on the scale of baryogenesis: Tbg > Aiss and Tew < Tbg < Aiss- 
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Figure 3. An example of a sequence of interactions present in the early universe which shifts the 
lepton number. 

• Tbg > Aiss 

(i) The Yukawa couplings of eq. (3.4) induce 2-body decays of the heaviest particle 
participating in each operator, e.g. N^ — t- Ei^ + ^*. Assuming these decays are 
kinematically allowed, they proceed with thermally averaged decay rate 

h]Mj (MA 

where Mj ~ Ajss is the mass of the decaying particle, hj is the relevant Yukawa 
coupling, and from now on Lorentz factors, such as [Mj/T), appearing in decay 



rates are understood to be present only for temperatures above the mass of the de- 
caying particle. To prevent washout one requires {nj/ni)T£) < i7 ~ g^''^T'^ /Mpi, 
for Tew < T < Tbg, where rij (rii) is the number density of species j (leptons). 
The strongest constraint arises from the latest time the density of the exotic par- 
ticles is still significant: at T ~ Ajss- This implies, 

,^ < (laV'^ ^1 '" - 10- (JLY'U^) '" , (3.12) 

■^ ~ V MpiJ VIOO/ VTeVy ' ^ ' 

where g are the relativistic degrees of freedom at T ~ Mj ~ Aiss . 
If this condition holds for yg, or h^ and h'c, or h^ and h' the lepton- number 
violating decays and inverse decays induced by these couplings are rare, and the 
BAU remains frozen. 

(ii) The thermally averaged rate of the decay mode X ^ Vl is, 

2 

Hence the lepton number of the universe is approximately conserved if: 



Mr?x^ 



V MpiJ VlOO/ VlTeVy ^ ^ 



(iii) The cross-section for the scattering process r/* + ,^* — t- <l> + $ is 



\2 1 



levr (r2 + Al 



(3.15) 



ISSJ 



and the corresponding scattering rate is F ~ n{T) a. In the relativistic regime, 
r oc T, and hence increases with temperature slower the Hubble parameter. In the 
T < Aiss regime, T becomes Boltzmann suppressed. Thus, for sufficiently small 
coupling 

this process is never in equilibrium. 

(iv) Lepton-number violating processes may also occur via some of the ISS degrees 
of freedom propagating off-shell. For example the scattering process Nji + x* ~^ 
Np —7- El + S, has an approximate cross section 

a-Mx)^^^^ (3 17) 

which leads to a bound 

<"^^^)^-'(ifo)'"(l^)"' 



Satisfying at least one of the bounds of eqs. (3.12), (3.14), (3.16) and (3.18) ensures no 
washout of the BAU, if baryogenesis has occured at temperatures Tbg > Ajss- However, 
as can been seen from eq. (3.9), these bounds would imply that the contribution to the 
active neutrino masses generated via this mechanism is insignificant. 

Tew < Tbg < Aiss 

The ISS degrees of freedom remain in thermal and chemical equilibrium even at tem- 
peratures much below the ISS scale, due to their gauge, Yukawa and scalar interactions. 
However, at Tbg < Aiss, their density is Boltzmann suppressed. This implies that the 
exotic leptons carry only a small fraction of the net lepton number of the universe. 

Consider for example again the 2-body decay rate of the ISS degrees of freedom via 
the Yukawa couplings of eq. (3.4). The decay rates at Tbg < Ajss are Tj « h'^Mj/lGn, 
where as before Mj is the mass of the decaying particle and hj is the Yukawa cou- 
pling that causes the decay. The lepton asymmetry will be approximately conserved if 
TjUj <C Hni, where n; ~ T^ is the lepton- number density, carried mostly by the light 
relativistic SM leptons, and Uj ss (MjT /Itt)^'"^ exp{—Mj /T) is the number density of 
the non-relativistic exotic leptons. 

For Yukawa couplings hj ~ 0{1), this yields the constraint 

Aiss 



Tbg 



> 42 , (3.19) 



(where we set Mj ~ Aiss)- The interaction considered in the case (ii) gives a similar 
constraint, while the scattering rates in cases (iii) and (iv) are sufficiently suppressed 
for: 

^ > 25. (3.20) 

Jbg 

We therefore uncover an important feature of this model: baryogenesis has to occur below 
the ISS scale. If it occurs before the EW phase transition, the ISS scale has to be sufficiently 
high, as seen by eq. (3.19). If it occurs after the EW phase transition, no limit on Ajss 
applies from considering the survival of the BAU today. Note that in this model, /x can only 
be generated after EW symmetry breaking (see eq. (3.8)), because the exotic fermions which 
gain a radiative Majorana mass carry hypercharge. This precludes resonant leptogenesis 
within this model, as leptogenesis must occur before the sphalerons switch off. It may, 
however, be possible that another of the ISS fields generates the asymmetry, with washout 
being suppressed by having the masses of at least one of the other exotic fields at a sufficiently 
high scale above Tbg- The construction of such a scenario, ensuring the satisfaction of the 
Sakharov conditions [56] and generation of the observed BAU, is beyond the scope of this 
work. 

3.3 Dark Matter 

In the Law/McDonald model the scalars, ^ and rj, and the fermions, N^_ and Nl, are odd 
under an accidental Z2 symmetry. The lightest of these is a DM candidate. Assume for now 
that the lightest state is in the scalar sector (we will briefly discuss the fermionic case below) . 
The neutral components of ^ and rj mix and the masses of the real and imaginary components 



split. The general mass squared matrix in the basis [Re{7]) Re{^^) Im{r]) /?7i(^'^)) is, 

Xi^,vl/A Ml 

V -A€<A,t;2/4 Ml j 

where we have taken \^(f,r] and [Xr]x to be real and positive without loss of generality as one 
can absorb the phase by redefining the fields, and 

Ml = 11%^ \r,^vl/2 + A^^i'2/2, (3.22) 

Ml = fil + Ag<^t;2 /2 + Xc)cvl/2. (3.23) 

The lightest of the eigenstates forms the dark matter. If the admixture of ^^ in the DM field 
is small (which is the case for M^ ^ M^) interactions with the Z boson will be suppressed. 
Scattering of DM with SM particles via the Z or Z' is also suppressed due to the mass 
splitting between the imaginary and real components of the DM field. (Gauge bosons interact 
off diagonally with the real and imaginary components of scalar fields.) Annihilation and 
scattering are dominated by the interactions with the Higgs fields, and for simplicity we will 
work in this regime. We then have a scalar Higgs portal DM candidate which has been 
studied extensively in the literature [57-60]. A possible difference from the standard Higgs 
portal arises when the heavier Higgs boson is light enough to influence the DM annihilation 
cross section [61-65]. We will explore this possibility in greater detail below in the context 
of the Law/McDonald model, though the discussion is largely applicable to generic Higgs 
portal models with an EW singlet Higgs mixing with the SM Higgs. 

3.3.1 Double Higgs portal DM 

We denote our (real scalar) DM candidate r], from now, and take its mass to be: 

where //^^^ > 0. One can choose parameters for the scalar potential (3.7) so $, and x gain the 
following VEVs: ($) = l/\/2(0 , v^Y , and (x) = %/\/^. The VEVs satisfy the fohowing 
equations: 

\^\^\\^x^l^\vl = ^. (3.26) 



To ensure that the potential is bounded from below and that (<^,x, ''?) = (^u)/v2,f^/v2, 0) 
is a minimum requires: 

A<^,Ax,A^>0, (3.27) 

\\^ < 4X^X^, (3.28) 

mIm > 0, (3.29) 

Ar?x ^ ~2y^A^A-^, (3.30) 

Aj7<^ > —2^/XrjX^. (3.31) 
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An additional condition for symmetry breaking is /x^ < 0, or //^ < 0, and to ensure there is no 
deeper minimum (at least at zero temperature) one can take fi'i > 0. The above conditions 
are necessary but not sufficient to ensure the potential is bounded from below. An additional 
sufficient condition is: 



1 J i^^<f>\x - '^^<)>x\ 



(3.36) 



The squared mass matrix for the physical Higgs bosons (j) and x is given by: 

One then finds the mass eigenstates {mH > ^n-/i), 

'^h = -^x4 + '^'t''"^ ~ V (-^x^x ~ ^'P'^iy + (^?ix^«'^x)^ (3-34) 

mjj = X^vl + X^vl + ^{Xx^l - A^v^ )2 + (X^^v^v^)^. (3.35) 

The mass eigenstates are superpositions of (j) and x, 

/i\ _ / cos^ sin^\ / 
^HJ \— sin 9 cos 6*^ yxy 

where the mixing angle 8 is given by the following, 

tan^ = hx^^^x (3 37) 

-^X^'x ~ ^-T^^^ + V (-^X^X ~ ^'P^D'^ + (A<i!>X^«'%)^ 

Note h (H) couples to SM particles the same way as the SM Higgs but with an additional 
factor of COS0 (— sin0) at the vertex. We identify h with the Higgs-like particle discovered 
at the LHC. 

In the limit where the H is very heavy the important DM interactions proceed through 
the SM-like Higgs, and one in effect returns to the standard Higgs portal scenario. If H 
is lighter, however, one now has annihilation to SM particles and interactions with quarks 
proceeding through both h and H propagators. We term this a double Higgs portal and we 
shall compare the two scenarios below. Similar ideas have been explored previously in the 
literature: refs. [61-63] consider two Higgs doublet models, while [64] considers fermionic 
DM with an exotic SM singlet Higgs. We consider scalar DM with an exotic singlet Higgs, 
a scenario which has been analysed previously in [65]. One important difference is that we 
focus more on the phenomenology of high DM masses, while the focus previously has been 
on the lower mass range Mjjm ^ 150 GeV. 

In terms of the Higgs mass eigenstates one finds the following interactions with the DM: 

—Cd -r] {Xfj^VujCosO + Xfj^v^sinO — idjj^sm0)h 

+ -T] {—Xfj^Viu sin 6 + Xr^^v^ cos 9 — ^r)x COS 6)H 

T] {Xqx COS 9 sin 9 — Xj^ff, COS sin 0)hH (3.38) 



r/2(A^0 sin^ 9 + A,,^ cos^ 9)H^ 



+ 7^^(A^0 cos^ 9 + Xr)x siii^ 9)h'^. 
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To reduce the number of parameters, and allow one to gain an appreciation of the phe- 
nomenology of the double Higgs portal, we study a special case where A,,;^; = A^,^ = A and 
f^vx ~ l-^lA'rjx' where we will vary A but keep /i^^ > and fixed. The interaction of the Higgs 
bosons with the DM then simplifies to: 

— £ D o ( ^'^ ^°^ ^ + ''^x ^^^ ^ ~ Sign[A] • /i„^ sin BjXr] h 



2 



1 



+ o ( "^x *^°^ ^ — Vw sin 9 — Sign[A] • /i^^ cos OjXrjH (3.39) 

+ \xr]^H^ + \xv^h^ 
4 4 

To see the effect of the second Higgs boson, we present two cases. Case A is where the 
heavy Higgs boson is taken to be heavy enough and the mixing negligible, so it is effectively 
decoupled. In this limit one obtains the standard Higgs portal. For case B, the double Higgs 
portal, we choose some parameters, namely Mh = 500 GeV, v^^ = 1 TeV, fi = 2 TeV, and 
examine the change in phenomenology for different choices of the Higgs mixing parameter 
sin 9. 

Furthermore, given the above choice of parameters, we have calculated A^, A,^, and X^f,^ 
as functions of sinO. One finds that condition (3.28) is satisfied for the entire range of the 
mixing angle (—71/2 < 9 < 7r/2), and that //?^ < so that symmetry breaking can proceed. 

As the heavy Higgs boson enters loop corrections to the W and Z boson propagators, 
a constraint on sm9 comes from the oblique parameters: T, S, and U. By following the 
analysis of [65-67] one finds that for Mh = 500 GeV the mixing angle is constrained to be 

|sin6'|<0.43 (3.40) 

by EW precision observations at 95% C.L. (assuming the other exotic particles give only 
negligible contributions.) 

3.3.2 DM relic abundance 

One can then calculate the annihilation cross section of DM into standard model particles. 
The relevant Feynman diagrams are shown in figure 4. To obtain the velocity averaged cross 
section, one should integrate over the velocity distribution. However a good approximation 
is obtained with, 

{aVrel) ^ icrVrel)\s=AMl^y (3-41) 

One can then simply compare to the required annihilation cross section, for the observed 
relic abundance [68], which for DM masses greater than 10 GeV is approximately [69], 

{cTVrel)relic ~ ^ >< 10"^^ Cm^/s, (3.42) 

where we have divided by 2 to take into account that we are dealing with only one degree of 
freedom for a real scalar rather than the two degrees of freedom for a Majorana WIMP. 

To obtain the desired M'^j^ with a choice /i^ > one wants to choose the required 
values of A to be negative (particularly for small DM masses), as can be seen from eq. (3.24). 
Some degree of fine tuning is of course required for small DM masses given our choice Vy^ = 1 
TeV. 

We have calculated the required coupling A = A^,^ = A,;;^, with the constraint A < 0. 
(The relevant expressions for the cross sections can be found in the appendix of ref. [60] 
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Figure 4. Possible DM annitiilation channels. The diagrams on the lower line give only a small 
correction and are not included in the calculations. 

for case A, and simple modifications for case B can be found by taking into account the 
Lagrangian in eq. (3.39).) The results for different choices of the mixing parameter are 
shown in figure 5. The observed relic abundance can be achieved with values of A ~ —0.1 
for most of the DM mass range. Note the falls in the required coupling at the points where 
2M£,M = Mfi/jj, as the intermediate Higgs boson is on-shell. 

We have also checked that the constraints of eqs. (3.30-3.32) can be satisfied with 
perturbative choices of the quartic coupling A,; (except for small regions of parameter space 
around sin0 ~ 0.05, where for low DM masses |A| > 1). 
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Figure 5. Required coupling for the correct DM abundance for different values of sin 6* and Case A 
(single Higgs portal) for comparison. Parameters have been chosen so that A = A^^ 



TeV, Mh = 500 GeV, and v^ = 1 TeV. 
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3.3.3 Direct detection 




Nuc Nuc 



Figure 6. DM Nucleon scattering through t-channel Higgs exchange. The effective coupUng of the 
Higgs with the Nucleon is indicated. 

DM particles will scatter off nuclei through t-channel Higgs boson exchange (see figure 
6). The spin independent DM- nucleon scattering cross section for our DM candidate r] is 
given by, 



A^{Mdm + Mnuc? \m]j 



ml 



where Mjsf^c — 0.94 GeV is the nucleon mass, fj\f parametrizes the effective Higgs-nucleon 
coupling, and the vertex factors are: 

A = -sin el ^ cos -sine - Sign[Al • 2^cos0|, (3.44) 

5 = COS 6l( COS 6* + ^ sin 6* -Sign [A] •2^sin6iV (3.45) 

Note that uncertainties in /jv result in roughly a factor of ~ 5 uncertainty in crf/j^j_j^^^ [70]. 
We use a value of fjy = 0.3 as indicated by lattice results [71]. 

Having calculated A, we then find f^^M-Nuc ^^'^ different values of sin 6*. The results 
are presented in figure 7. XENONIOO exclusion limits from 255 live days of data taking 
[72], labeled XENONIOO, and projected XENONIT exclusions for 2017 in the case of no 
WIMP-nucleon scattering [73], labeled XENONIT, are plotted for comparison. We see that 
Mdm ^ 80 GeV is still largely viable, but that much of the remaining parameter space can 
be ruled out by XENONIT. 

Note that due to different interferences the scattering rate can be above or below the 
zero mixing case. A minimum is reached for sin0 = —0.052, which is exactly where, 

A B\ , , 

+ = 0, (3.46) 

rriu mj^J 

so that the DM-nucleon cross section is zero. Similar interference was noted in the case of a 
two Higgs doublet extension of the Higgs portal in ref . [62] . 



14 



10- 



10' 



b 10- 



10- 



10- 





I 








41 
43 
4^ 




it^ '"Sj. 


^ 


-i 




^ 


/^ — z 


— 




/-- — 




47 


■ 


-- w 




■ 


49 




T 







10 



20 



50 100 

Mdm [GeV] 



200 



500 



1000 



Case A 

sin (6) = 0.1 

sin (6») = 0.15 

sin (6) = -0.1 

sin(0) = -O.15 

XENONIOO 

XENONIT 



Figure 7. DM- Nucleoli scattering cross section for case A and case B with different values of sin(6'), 
the upper limit from XENONIOO, and the projected sensitivity for XENONIT in 2017. 

3.3.4 Higgs signals at colliders 

The double Higgs portal can be tested in a number of ways at colliders [62, 63, 74, 75], some 
of them common with the usual single Higgs portal [67, 70, 76, 77]. Mixing with the exotic 
Higgs will reduce the production of the SM-like Higgs, 



a. 



pp—>h 



COS 



SM 



(3.47) 



The total width will also be reduced by a factor of cos^0, however the SM-width for rrih = 125 
GeV, roughly 4 MeV, is already far below current experimental resolution. The discovery of 
a SM-like Higgs boson at the LHC [78, 79] allows us to constrain this scenario. The observed 
signal strength over the SM expectation and la uncertainty is [80, 81]: 



o-/crsM = 0.88 ± 0.21 (CMS), 

a/asM = 1-3 ± 0.13 (stat.) ± 0.14 (syst.) 



(ATLAS). 



(3.48) 
(3.49) 



Combinations of the results put a limit on the universal suppression of all channels so that: 
(j/(JSM ^ 0.8 at the 2a level [82-84]. This translates into a limit. 



sin 6*1 < 0.45, 



(3.50) 



if one assumes no invisible decays of the 125 GeV resonance. This is comparable to the 
constraint from EW precision observables, eq. (3.40), however the limit on the mixing angle 
from measurement of the signal strength is generally applicable and not limited to our choice 
Mh = 500 GeV. 
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If the DM is light enough, the Higgs can decay invisibly to two DM particles. The 
invisible partial decay width of the Higgs {h — t- rjr]) is given by, 



A^ / . , ' \2 / 4m^ 



1/2 

r = — f Vuj cos + v^^ sin 6 — Sign[A] • 2^^^ sin 9] I 1 ^ 1 . (3.51) 

This leads to a large invisible branching fraction Br{h — )• ryr]) > 0.99 for light DM masses 
Mdm ^ 10 GeV, except for a very narrow region around a mixing angle sin^ = —0.049 
where the effective coupling of the hT^rj term is zero. This helps to rule out such light DM 
candidates as the direct detection experiments lose sensitivity for low mass DM candidates. 
If the mixing angle is large enough, the heavy state could also be detected. For a Higgs 
mass of 500 GeV, CMS data currently sets a 95% C.L. [85]: 

a/asM ^ 0.21. (3.52) 

The most stringent limit on the mixing angle will come from assuming 2Mdm > Mh, so 
there is no invisible width. There is also an additional complication oi H ^ hh decays [86] , 
however, for the parameters we have chosen the relevant branching ratio turns out to be 
negligible. The limit on the mixing angle from high mass Higgs searches is then: 

I sin 6*1 < 0.46. (3.53) 

Thus the most stringent limit on the mixing angle in case B with Mh = 500 GeV comes from 
considering EW precision observables, eq. (3.40), while the most widely applicable constraint 
comes from measurements of the signal strength of the 125 GeV resonance, eq. (3.50), as no 
assumptions about the mass of the heavier Higgs or VEV have to be taken as an input. 

3.3.5 Fermionic DM 

Alternatively one may consider fermionic Higgs portal DM. The important terms are then: 

- £ D h^xN^iNRf + h'^xNliNLT + MnNINr + H.c, (3.54) 

where the lightest mass eigenstate of the A'^^;,, {NrY admixture is the DM candidate (pro- 
tected by the same accidental Z2 as in the scalar DM case). In such a scenario DM is 
constrained to be very heavy from direct detection experiments, M^m ^ 2 TeV, unless the 
DM mass is at one of the two Higgs resonances, there is parity violation, or a Sommerfeld 
enhancement [64]. So while we have studied the scalar DM case above, either scalar or 
fermionic DM is viable in the Law/McDonald model. Further generalizations taking into 
account interactions with the Z' in the fermionic DM scenario are also possible, but we will 
not pursue the details here. 

3.4 Concluding remarks on the Law/McDonald model 

We have seen that in the Law/McDonald model, baryogenesis is severely constrained. The 
explicit lepton number violation will tend to wash out any baryon asymmetry created at 
a high scale. However at temperatures roughly a factor of 25-40 below the ISS mass scale, 
Boltzmann suppression of the thermally averaged reaction rates means that the ISS fields will 
not contribute to washout (though not necessarily all of the ISS fields must have such high 
masses) . Unfortunately the small L violating parameter is only generated after EW symmetry 
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breaking, meaning that resonant leptogenesis is not possible in the Law/McDonald radiative 
ISS. 

We found that the Law/McDonald model can easily accommodate the observed DM 
abundance through the Higgs portal. Also, in a simple extension, an additional Higgs state 
can also contribute to the DM annihilation rate and DM nucleon scattering rate. Apart from 
narrow regions of parameter space, either on resonance points for the annihilation rate, or 
near complete destructive interference for the DM nucleon scattering rates, these models can 
for the most part be tested by the XENONIT experiment. 

4 Ma radiative inverse seesaw model 

4.1 Review of the model 

We now consider the Ma radiative ISS model [29]. An important difference in this model is 
that {B — L) is a protected symmetry down to the ISS scale. First let us review the main 
features of this model. It employs an additional U{1) which originates from a grand unified 
theory (GUT) with the following breaking pattern, 

SO{10)^SU{5)C^U{i}^ 

-^ SU{^)c SU{2)l C/(l)y f/(l)x. ^ ' ^ 

The exotic charge is given by the following, 

Q^ = 5{B -L)-AY = b{B - L) + AI^l - ^Qem- (4.2) 

There are two scalars, r/i and ry2, that transform under U{1)^^ as per 

Qx(m) = l, QxM = 2, (4.3) 

and are neutral with respect to the SM gauge group. These will be used for the radiative 
ISS and to break the new C/(l)^. To allow the usual Yukawa interactions, 

C = XdQl ■ ^dR + Xue^'Ql^^luR + Xeh • '^cr + Xue'^'lTa^lNR + H.c, (4.4) 

the SM Higgs doublet must carry a U{1)^ charge, 

$ = (<A+,0°) -(1,2,1/2, -2). (4.5) 

The fermion multiplets of 50(10) contain a right handed neutrino, Nr ~ (1, 1,0, —5), per 
generation of SM fermions. To these are added more SM neutral fermions with Qy. given by, 

lxS3L~-3, 4x52L~2, 5xSiL^-l, (4.6) 

where the required number of copies has been indicated to form an anomaly- free set. The 
following Yukawa terms involving the exotic scalars are permitted, 

C D fRsN^S^Lvl + f23'(S^'S2LVl + fi2(S^SiLvl + /ii(S^5il7?2 + H.C. (4.7) 

The scalar potential is given by, 

V = {Ar]l4 + H.c.) + Xnivlmf + A22(r?^r/2)^ + Xi2ir]lm)i4m) 

+ /i|$"f$ + A$($t$)2 + Ai$($"f$)(7?|r?i) + X2<s>{<^^<^){vlm) (4.8) 

+ f^iivlvi) + fj-livlm), 
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where ^ is a coupling constant with dimension of mass. The U{l)y. symmetry is spontaneously 
broken by, 

im) ^ vjV2 / 0, im) = 0. (4.9) 

Before radiative corrections are taken into account the neutral fermion mass matrix is given 
by, 

/ niD \ / {ulT \ 



^=7;{W {NrY Ssl S,l) 



mo Mn 
Mn 
V mJ 



Nr 

iS3Lr 



(4.10) 



where rriD = ^i^v^/V^, Mn = /fis^/v^, and Mi = 1/2/11%- Note the S2L are massless at 
tree level, and that the scalar interaction term, 



Arjirjiry^ + A*ri['q[ri2 — > Arjiriiv^^/\/2 + yl*r?|r/ji;^/\/2, (4.11) 

induces a mass splitting between the mass of the real (m/j) and mass of the imaginary (m/) 
components of r/i , 

m\-m] = 2V2Av^, (4.12) 

where we set A to be real and positive without loss of generality as any complex phase can 
be absorbed into the fields. 
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Figure 8. Radiative mass generation for 5*2 and 53 in the Ma model. 

The exchange of Re{rji) and Im{rii) then give one-loop radiative Majorana masses to 
the S2L as can be seen in figure 8: 



Mo 



(/l2)'Mi 

167r2 



mr 



m 



R 



-^lnf4 
Ml \Ml 



m. 



rln 



TJIt 



Ml \Ml 



nil 



(4.13) 



Siinilarly one obtains two-loop radiative Majorana masses for the /Ssl, 



f^^Ms 



(/23)'M2 

167r2 



m 



R 



mj^- M2 



2 v^l 



rri' 



r^lnf4 
'-M| \Ml 



(4.14) 



which are exactly the small parameters required in the mass matrix of eq. (4.10) for the ISS. 

Note rji and S2L are odd under an accidental Z2 symmetry, while all the other particles 

are even. We will see below, that for typical choices of the parameters, S2L is lighter than rji 
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as its mass arises radiatively and therefore S2L is the DM candidate. One can find various 
approximate forms for the radiative masses depending on the sizes of various parameters. If 



we take mj^ S> mj^ 



m. 



2\^Av^,Mi, eq. (4.13) can be approximated as, 



M2 



V2{fi2?Av^Mi 



87r2 



m, 



(4.15) 



while for natural choices of the parameters m'j^,m'j ^ M2, one obtains from eq. (4.14): 



^^ 



M, 



ifl2)Hf23?AXMl 



327r4 



rrij 



(4.16) 



To give a numerical example, for v-^ ~ Mi ~ A ~ 0(10 TeV), and mj = 50 TeV, one finds 
M2 ~ (/i2)^ X 7 GeV and M3 ^ (/23)^(/i2)^ x 5 MeV (in very close agreement to what the 
exact formulas give). 

4.2 Constraints from BAU washout 

Until f^(l)x is broken, (B — L) is an exact symmetry of the Ma model, and generation or 
washout of the BAU will not take place. ^ The VEV v-^^ breaks Q-^ = b{B — L) — 4Y, and 
hence the (B—L) as defined above. Baryogenesis can then take place, either through resonant 
leptogenesis, or some other suitable baryogenesis scenario with Tbg ^ T-^ ^ v^. However, as 
shown in figure 9, the ISS fields can now also wash out the asymmetry. 







<m> 
I 

t 



m 




Figure 9. Washout process involving the exotic fermions and scalars. Also shown are the Yukawa 
interactions which link this washout process to the SM sector. Exchange of the real and imaginary 
components of 771 leads to destructive interference. Not shown are the interactions involving /12 and 
/ii which ensure S2L cannot be assigned non-zero SM lepton number. 



A new conserved global {B 
following parameters vanishes: 



L) can be defined, even with (7/2) 7^ 0, provided any of the 



^u, fm, f23, fl2, /ll, A. 



(4.17) 



^This constraint could be circumvented through sequestration of positive and negative B — L charge into 
the visible sector and a hidden sector. This is a popular mechanism in asymmetric DM models, which aim 
to dynamically relate the visible and dark matter relic abundances [87-99]. In the context of the model 
studied here, it would require introducing additional fields and interactions, such that the B — L symmetry 
encompasses two low-energy global U(l) symmetries - the {B — L)v carried by the SM fields and a dark baryon 
number Bd- The construction of such a model in an S'O(IO) GUT is beyond the scope of this paper. 
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If the interactions mediated by one of the above parameters are out of equihbrium the BAU 
is safe. We first examine what values the parameters must take so as to not wash out the 
BAU. In section (4.3) we examine the possibihty for resonant leptogenesis, where the hghtest 
pseudo-Dirac pair generates the BAU. 

• Tbg » ^iss ~ % 

(i) The decays and inverse decays Nr o //, + $ are controlled by X^ . The decay rate 
is approximately, 

r.^(^). ,4.18) 

requiring nj\fT < riiH where n^ (ni) is the Nji (lepton) number density leads to 
a washout avoidance condition: 

(ii) The mixing between the Nji and ^3^ states is controlled by fji^ = ^plM^ jv^^. 
Scattering $ + /^ — )• Nji/ S'n -^ Vi'^ ^21 has an approximate cross section, 

^ ^ (A./.3)^ Ml 
The scattering rate F ~ aT^ is inefficient compared to H provided: 

and for such values washout is avoided, 
(iii) The Yukawa couplings /23 and /12 induce two body decays of the heaviest par- 
ticle involved in the relevant interaction (such decays are kinematically allowed 
given S2L-, which has a radiative and hence negligible mass compared to the other 
particles, can appear in the final state). We estimate the decay rate as. 



167r V T 

where fi denotes the relevant Yukawa coupling and Mj the mass of the heaviest 
particle. This leads to a washout avoidance condition of: 

(iv) Scattering processes S2L + '?i — ^ Si/ Si — )• S2L + "Hi can be suppressed with small 
/ii = M\l \f2v^. We estimate the cross section to be: 

[h2f Ml 
This gives a washout avoidance condition: 
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(v) Washout may proceed through decays and inverse decays t]i provided they are 
kinematically allowed. However, for a small enough A, and hence small enough 
mass splitting between Re{r]i) and Im{r]i), destructive interference leads to washout 
suppression. We estimate the decay rate as, 

r^ (|/23p + |/l2p)M, M,^ 

16ir T 

and the effective washout parameter to be, 

• (5^ (4.27) 



A'., :l 



T=Mi 



where 6 = {m^ — mi)/T. One finds this translates into a washout avoidance 
condition: 

A^^^-'O-^-AS^^W'ijI^f^^)"^. (4.28) 

If such decays are not kinematically allowed washout can instead proceed through 
a scattering process mediated by off-shell r/i. We estimate the cross section as: 



^ ^ ...o.^...u , __^ _ __ .\. ^4_29) 



(|/2 3| + |/l2|)V 1 ^^ ^2 

1287r \T'^ + m\ T"^ + m] , 

(I/23I + I/12I)" / 2V2An 



2 
^ * T^. (4.30) 



1287r I (r2+m|)(r2+m2) 

This translates into a washout avoidance condition: 

Satisfying any of the above washout bounds means the BAU is safe. However, simi- 
larly to the Law/McDonald model, such small choices for the couplings ruins the ISS 
mechanism. 

Tbg < Aiss 

So as to ensure the preservation of the BAU, while allowing for the ISS mechanism, 
we may suppress washout by taking one or more of the exotic particle masses to be 
higher than Tbg- This ensures that at Tbg the number density of that particle will 
be sufficiently Boltzmann suppressed for it to not play a role in washout. For decay 
processes the requirement on the mass scale is: 

^>42, (4.32) 

while for scattering processes the constraint is relaxed to roughly, 

Aiss 



Tbg 



> 25. (4.33) 



This is again similar to the washout in the Law/McDonald model. For Tbg < Tew the 
bounds of eqs. (4.32) and (4.33) do not apply as the sphalerons are no longer active. 
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4.3 Resonant leptogenesis 

For typical parameter choices, the hghtest pseudo-Dirac pair of neutrinos, Na wih generate 
the BAU in resonant leptogenesis. There is an interplay between the CP violation and 
washout, which both depend on fi, the mass splitting of the Na pair. So as to show that this 
model can generate the right order of magnitude for the parameters, let us take as a guide the 
example of [28], which found the observed BAU can be generated with 6 = fi/Ta ~ 10~^ (Fq, 
is the width of the N^ and fi the mass splitting of the pair), heavy neutrino mass Mjsfa ~ 0(1) 
TeV, and Yukawa coupling Aj^ « 5 x 10~^. 

The scalar mass mj can be taken to a sufficiently high scale so that the other interactions 
of the ISS do not interfere with the generation of the BAU. Let us take mi ~ 50 TeV and 
v^r^ Ar^ Ml r^ W TeV, then Mn ~ 0(1) TeV can be achieved with /rs ~ 0.1. Note that 
solutions of the Boltzmann equations taking into account decay and scattering processes show 
the BAU is generated at Tbg ~ -^Afa/10, which allows the out-of-equilibrium condition for 
Na to be met even with a relatively low Z' mass: Mz' ^ 2Mjvq [54]. 

Choices of 6 = 10-^ A^ = 0.05, and M^a = 2 TeV imply: /i = 0.2 x 10"^ GeV. This is 
compatible with a radiative origin in the Ma model as, for these parameters, eq. (4.16) gives 
fJ- ~ (/i2)^(/23)^ X 5 X 10"^ GeV. The light neutrino mass is rriu ~ 10~^ eV. 

In the (5 <C 1 regime, the CP asymmetry of N^ decays is [28], 



A^ 



fin 



IGtt {fi^ + fi^ + 2fifi cos a 



sma. 



(4.34) 



where a is a physical phase of the ISS mass matrix and fi is the Majorana mass fiNii{NiiY. 
In this model one typically expects fl ~ 10"'^ /i, due to suppression by additional Yukawa 
couplings in the radiative generation of /i, which allows for e ~ 10~®. The washout factor is 
-fi'eff ~ 6. Taking into account a sphaleron reprocessing and dilution factor, the final baryon 
number density to photon number density is, rfs ^ 10~^|e|, and can easily be made to match 
the observed value tjb ^ 6 x 10"^*^ [100] in this scenario. 

A more detailed study finding the full form of the 3x3 matrices rriD and fi, and mass 
spectrum for the pseudo-Dirac pairs compatible with the observed oscillation parameters, 
BAU, and 50(10) origin of the model is beyond the scope of this paper. 



4.4 The p parameter 

In the Ma model, the Higgs doublet responsible for EW symmetry breaking also carries U{l)y- 
charge. This leads to a tree level correction to the Z boson mass, and consequently to the p 
parameter. 



where the Weinberg angle is defined in terms of the weak isospin, g, and weak hypercharge, 
qy, coupling constants as tan9]y = gy/g- Experimental observations give a value p — 1 = 
4^*^^ X 10~^ [101]. The masses of the neutral gauge bosons are given by. 
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g'^+gl + IQgl 



1 + 
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V 



+ 4[45x5'y sin Ow + ^9x9 cos OwV 



(4.36) 
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so to remain within the 2a hniit of the p parameter, one requires v^ > 7.7 TeV (given 
9x ~ 9y)" The Z — Z' mixing angle is given by, 

, „, 2{g^gY sin Ow+ggx cos Ow) , . „„s 

SO one finds |</>z| — (^«>/%)^ ^ 10~^ given u^ > 7.7 TeV. We will see the consequences of 
this for our DM candidate below. 

So far we have dealt with mixing coming from the mass matrix. As the theory contains 
particles charged under both C/(l)y and C/(l)^, kinetic mixing between the respective field 
strength tensors will be generated radiatively. The full mass matrix for the Z and Z' bosons, 
taking into account both kinetic and mass mixing, can then be found [105]. By following this 
procedure we have checked that including such a term makes no qualitative difference to our 
discussion above. Furthermore one typically expects radiatively generated kinetic mixing to 
be small [106] and therefore negligible for mixing in this model. 

There is also a competing limit from the LHC. For example, using the ATLAS limit on 
dilepton resonances [107] and a choice of g-^ = 0.1, one finds a lower bound of M^/ > 2.7 
TeV. Given that Mzi — 2g^v^, this corresponds to \(j)z\ ;$ 6 x 10^^. Larger g^ means 
stronger collider limits on Mz' ■ However, ATLAS and CMS currently only provide limits for 
Z' resonances up to 3.5 TeV [107, 108], and as the limit on the possible signal strength is 
rising for high mass regions, we do not pursue the details here. Suffice to say the bound on 
the mixing angle is approximately [(/'zl i$ 10^^. 

4.5 Dark matter 

4.5.1 DM candidate mass 

In the Ma radiative ISS, there exists an accidental Z2 symmetry under which only rji and S2L 
transform. The lightest of these species is stable and forms a DM candidate. By examining 
the scalar potential in eq. (4.8) we see that the masses of the real and imaginary components 
of rji after symmetry breaking are given by, 

ml = (/i^i + Ai0t;2 /2 + \uvl/2f + V2Av^ (4.38) 

mj = inl^ + Ai0t;2 /2 + Ai2z;^/2)2 - V2Av^, (4.39) 

where A > and for reasons of naturalness one typically expects /x^i ~ f -^ ~ A. The fermion 
S2L gains a radiative Majorana mass, M2, at 1-loop order, and one typically expects, 

m/ > M2, (4.40) 

if fine-tuning is absent. To provide a sense of the typical values for M2, we have have plotted 
M2 as a function of ^VAv^, with Mi = ^/Av-^ for different choices of m[ in figure 10. One 
can see that as one moves away from the fine tuned limit mj = 0, the radiative mass M2 
decreases to be in the range, 

M2 « (10-2 - W-^){h2)\. (4.41) 



■^The constraint can also be satisfied with g-^ < Qy /95 and v-^ < 5 TeV. Given the GUT origin of the 
model, we consider such small choices for g^ unattractive. Furthermore the interactions of the resultant light 
Z' are constrained to be so weak by observations of Caesium transitions [102-104], that they would also lead 
to overclosure of the universe. 
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Figure 10. A/2/(/i2)^ as a function of ^,/Av^, with Mi = ^JAv^, for different choices of m/. 

4.5.2 Cold DM scenario 

Given the discussion of the p parameter one requires v^^ Vy^- We will discuss the DM relic 
abundance in light of this constraint here. Let us assume the S2L species have a thermal 
abundance in the early universe, as will be the case unless g^^ < 10~^ <C gy- There is 
furthermore an absence of any interactions of the form 5'2l5'2l(o") where a represents any of 
the scalar fields in the theory. As S2L is charged only under C/(l)^, it can annihilate only 
through the Z' , and due to mixing of the gauge bosons, through the Z boson [109] (see figure 
11). 




/; 



SM 



/; 



SM 



Figure 11. Annihilation process of S2L which determines its relic density. 



So as not to overdose the universe one requires the annihilation cross section to be 
larger than {(t\v\)dm ^ 1-9 x 10^^ GeV^ [69]. However annihilation through the Z boson is 
suppressed by the small Z — Z' mixing angle. For example, consider the 2M2 < Mz regime. 
It suffices for our purposes to use a simple estimate for the annihilation cross section: 



{a\v\) 



sin^ (l)zGlMl 

2^ 



(4.42) 
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where (j)z is the Z — Z' mixing angle and Gp is the Fermi constant. Requiring (o"|w|) > 
{(^\v\)dm, ^-iid substituting for the mixing angle 4>z — (vw/v^)'^, one finds: 

Gf \Vw/ 

Given the constraint from the p parameter, v^> 7 TeV, this demands M2 > 7.5 TeV, which 
contradicts our original assumption of 2M2 < Mz- One also finds a too small cross section 
in the 2M2 > Mz regime. 

Annihilations occuring through the Z' are suppressed by the large mass M^/. So anni- 
hilations through the Z' are also too weak. 

Finally we note that the cross section will be much higher at a resonance point, where 
2M2 — Mz or 2M2 — Mz' [109], and this could avoid overclosure. However, this corresponds 
to a fine tuning of the DM mass. 

From this discussion we see that for a simple cold DM scenario the S2L DM candidate 
overdoses the universe if we demand no fine tuning. One would have to go to more convoluted 
DM scenarios such as warm dark matter to perhaps find viable regions of parameter space. 
The simple cold DM case is ruled out. 

4.6 Concluding remarks on the Ma model 

The Ma radiative ISS model has several positive features for baryogenesis and washout avoid- 
ance. Resonant leptogenesis may be accommodated, as the mass splitting of the heavy neu- 
trinos is generated before the EW phase transition, allowing the interplay between heavy 
neutrino decays and sphalerons crucial for leptogenesis. Also, washout due to the other AL 
interactions required for the radiative ISS can be suppressed with suitable choices of mass 
parameters. 

On the other hand we have also seen a drawback of the Ma model with respect to 
cold DM. We found that the constraints from measurements of the p parameter translated 
into DM annihilations through the Z and Z' gauge bosons which are too weak to avoid 
overclosure. These simple cold DM scenarios are therefore ruled out for the Ma radiative ISS 
(apart for possibly small areas of parameter space corresponding to resonant enhancements). 
One would therefore have to go to more convoluted scenarios such as warm dark matter to 
find a suitable DM scenario. 

5 Conclusion 

Motivated originally by how new physics explanations for the baryon asymmetry, dark matter, 
and neutrino oscillations can overlap with each other, we began by looking at the properties 
of the inverse seesaw. We saw how the inverse seesaw mechanism offers various advantages, 
namely a scale low enough for direct experimental tests, and the possibility of resonant 
leptogenesis, or at least a suppressed washout rate if some other mechanism is responsible 
for the baryon asymmetry. 

The inverse seesaw mechanism does comes with a cost: the small lepton number violating 
parameter that must enter the mass matrix. Despite its technical naturalness, there may be 
a more appealing mechanism to generate it. One of the options is a radiative origin. We 
therefore continued our investigation by focusing on two radiative inverse seesaw models. 

It was found in the case of the Law/McDonald radiative inverse seesaw model, the DM 
abundance is easily explained with the Higgs portal mechanism, or the double Higgs portal. 



25 



The advantages of the inverse seesaw with respect to baryogenesis, however, are lost due to 
the extra fields present which break lepton number exphcitly with their interactions. In case 
of the Ma radiative inverse seesaw model, the situation is reversed: it nicely accommodates 
resonant leptogenesis, but does not feature a simple cold dark matter candidate. 
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